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The elusive Heisenberg limit in quantum- 
enhanced metrology 
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Quantum precision enhancement is of fundamental importance for the development of 
advanced metrological optical experiments, such as gravitational wave detection and 
frequency calibration with atomic clocks. Precision in these experiments is strongly limited by 
the l/V/V shot noise factor with N being the number of probes (photons, atoms) employed in 
the experiment. Quantum theory provides tools to overcome the bound by using entangled 
probes. In an idealized scenario this gives rise to the Heisenberg scaling of precision 1//V. 
Here we show that when decoherence is taken into account, the maximal possible quantum 
enhancement in the asymptotic limit of infinite N amounts generically to a constant factor 
rather than quadratic improvement. We provide efficient and intuitive tools for deriving the 
bounds based on the geometry of quantum channels and semi-definite programming. We apply 
these tools to derive bounds for models of decoherence relevant for metrological applications 
including: depolarization, dephasing, spontaneous emission and photon loss. 
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uantum- enhanced metrology aims to exploit quantum fea- 
| tures of atoms and light such as entanglement, for meas- 
uring physical quantities with precision going beyond 
the classical limit 1-4 . A prominent example is that of an optical 
interferometer, where interference of photons at the output port 
carries information on the relative optical path difference between 
the interferometer arms. When standard laser light is used, the 
observed results are compatible with the claim that each photon 
interferes only with itself' 5 and the whole process may be regarded 
as sensing with N independent probes (photons^Parameter estima- 
tion with N independent probes yields the l/y/N standard scaling' 
(SS) of precision 6 . Entangling the probes however, can in princi- 
ple offer a quadratic enhancement in precision, that is, the l/N or 
'Heisenberg scaling' (HS) 7-11 . Such strategies have been experi- 
mentally realized in optical interferometry 12-16 with exciting appli- 
cations in the quest for the first direct detection of gravitational 
waves 17 ' 18 . Moreover, the same quantum enhancement principle 
can be utilized in atomic spectroscopy 19,20 where the spin-squeezed 
states have been employed for improving frequency calibration 
precision 21-24 . Alternative approaches to beat the SS without 
resorting to quantum entanglement include multiple-pass 25 ' 26 and 
non-linear metrology 27 ' 28 . 

Unfortunately, both the theory 29-35 and experiments 36 ' 37 con- 
firmed the fragility of the above schemes when noise sources such 
as decoherence are considered, and it has been rigorously shown 
for particular models that asymptotically with respect to N, even 
infinitesimally small noise turns HS into SS, so that the quantum 
gain amounts to a constant factor improvement 38-40 . 

In this paper, we develop two methods that allow us to extend 
these partial results to a broad class of decoherence models, and in 
particular to obtain fundamental bounds on quantum enhancement 
for the most relevant models encountered in the quantum metrol- 
ogy literature. The two approaches complement each other in terms 
of provided intuition and power. The first method elaborates on 
the idea of classical simulation (CS) 41 and provides a bound based 
solely on the geometry of the space of quantum channels (SQC). 
When applicable, it gives an excellent intuition as to why the HS is 
lost in the presence of decoherence, but fails to yield useful bounds 
for some relevant decoherence models. The second method is based 
on the channel extension (CE) strategy 42 and requires the optimi- 
zation over different Kraus representation of a quantum channel. 
However, unlike the earlier work 40 , which involved making an edu- 
cated guess about the appropriate class of Kraus representations, 
our bound can be cast into an explicit semi-definite optimization 
problem, which is easy to solve even for complex models. The 
power of the method is demonstrated by obtaining new bounds for 
depolarization and spontaneous emission models, and re-deriving 
asymptotic bounds for dephasing and lossy interferometer with 
unprecedented simplicity. 
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Figure 1 1 Quantum metrology and the CS idea, (a) General scheme for 
quantum-enhanced metrology. N-probe quantum state fed into N parallel 
channels is sensing an unknown channel parameter (p. An estimator <p is 
inferred from a measurement result on the output state, (b) CS of a quantum 
channel. The channel is interpreted as a mixture of other channels Ax, 
where the dependence on cp is moved into the mixing probabilities p<p(X). 



We pursue the estimation problem by restricting our attention 
to estimators, which are unbiased in the neighbourhood of some 
fixed parameter value, for which the quantum Cramer-Rao bound 
(CRB) holds 



(2) 



where Fq is the quantum Fisher information (QFI) 9 . Maximiza- 
tion of QFI over input states sets the limit on the achievable 
quantum-enhanced precision: 



F N = max F f 



~N 



(3) 



The states that maximize the QFI and yield the HS for decoherence- 
free unitary channels are typically highly entangled: the GHZ state 
in the case of atomic spectroscopy 43 , and the NOON state in the case 
of optical interferometry 44 . In the presence of decoherence, the opti- 
mal input states do not have an intuitive form and the maximization 
of the QFI (with rising N) becomes hard even numerically 20 ' 29 ' 30 . 

The typical behaviour of the estimation uncertainty in the 
presence of decoherence is depicted in the log-log scale in Fig. 2, 
showing that asymptotically in N the quantum gain amounts to 
a constant factor improvement over the standard 1/ViV scaling 
achievable with independent probes. The key result of this paper 
is to provide a method for a general and simple calculation of this 
constant factor improvement. 



Results 

Bounds on precision in quantum-enhanced metrology. The typical 
quantum metrology scenario is depicted in Fig. la. An ensemble of 
N quantum systems undergo in parallel the same transformation 
A<p, which depends on an unknown physical parameter (p. The 
output state is measured and the outcome is used to compute an 
estimate (p of the parameter <p as summarized below 

„->A™[y]->*. (i) 

The task is to find the optimal (possibly highly entangled) input 
state p N and the most effective measurement strategy to minimize 
the estimation error A(p N . Note that since the decoherence process 
is assumed to act independently on each of the probes, the global 
channel is described by the tensor product A® N . 

2 



Classical simulation. To understand the idea of CS, we need to 
think of quantum channels in a geometrical way 45 . A quantum 
channel is a completely positive, trace-preserving map acting on 
density matrices. The space of all such transformations is convex: if 
A, A' are two channels, then pA + (1 -p)A' can be realised by ran- 
domly applying A or A' with probabilities p and 1 -p. The channels 
that cannot be decomposed into a convex combination of differ- 
ent channels (for example, the unitary transformations) are called 
extremal'. Note that while all interior points of the SQC are non- 
extremal', the boundary contains both the extremal as well as some 
non- extremal channels. 

We say that the family A^ is classically simulated' 41 if each 
channel is written as a classical mixture of the form 

A (p [p] = jdx p (p (x)A x [p], (4) 
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Figure 2 | Estimation precision in presence of decoherence. Log-log plot 
of a generic dependence of quantum-enhanced parameter estimation 
uncertainty in the presence of decoherence as a function of the number of 
probes used. While for small number of probes the curve for achievable 
precision follows the Heisenberg scaling, it asymptotically flattens to 
approach const/ V/V dependence. The 'const' represents the quantum 
enhancement factor. The exemplary data corresponds to the case of phase 
estimation using N photons in a Mach-Zehnder interferometer with 5% 
losses in both arms. 



where the unknown parameter enters only through the probability 
distribution p™ of a random variable X that indicates which channel 
to pick from the set {A x } (Fig. lb). If X h . . .,X N are the independent 
hidden random variables used to generate the parallel channels, we 
can rewrite the estimation problem as 



>q>. 



(5) 



Since conditionally on the values of Xj the output state does not carry 
any information about <p, the estimation precision in the above sce- 
nario is at most equal to that of the classical problem of estimating (p 
given N independent samples from the distribution p^ 



(6) 



Hence, using the classical CRB 46 , we obtain a lower bound on the 
uncertainty of the original problem equation (1) 



-\dx 



PcpM 



(7) 



where F c \ is the classical Fisher information; see Methods for an 
alternative derivation. 

Provided that py{x) satisfies regularity conditions required in the 
derivation of CRB 46 , the SS of precision follows immediately. More- 
over, any particular CS yields an explicit SS bound on precision. If 
many CSs are possible, we obtain the tightest SS bound for A(p N by 
choosing the worst' decomposition yielding minimal F c \[p^[ , as shown 
below. On the other hand, HS is possible only when the above condi- 
tions are not satisfied. This happens, for example, in the decoherence- 
free case where are unitary channels, which are extremal points 
of the SQC and the only admissible p<p in equation (4) is the irregular 
Dirac delta distribution being zero on all channels except A^. 

The QFI at a given <p 0 depends only on the output state and its 
first derivative at <p 0 . This implies that any family of channels A^, 




Figure 3 | Local classical stimulation. Schematic representation of a local 
classical stimulation (CS) of a channel A<p at q>Q that lies inside the convex 
set of quantum channels (solid oval). The optimal CS has to be valid only 
in the neighbourhood of A~ along the curve f A m \ (solid arched line) 

and corresponds to a binary mixture of channels A+, which rest on the 
tangent (dashed) line at the two outermost points. Then, the precision of 
estimation can be lower bounded just using the distances e+. 



which 'locally coincides' with the original one, that is, 
[p] = [p] , 3^ [p]|^ o = 3^ [p]| 



(8) 



achieves the same maximum QFI in equation (3). It is therefore 
enough to consider the 'local classical simulations', that is, any mix- 
tures reproducing the channel and its first derivative at given (p 0 . 
As proven in Supplementary Methods, the CS with the smallest F c \ 
can be constructed using two channels, {A + , A_}, which lie on the 
tangent line to the 'channel trajectory' at the two outermost points 
situated on the boundary of the set of channels (Fig. 3). The local CS 
around <p 0 reads explicitly 



Ap[p] = pjA + [p] + PpA_[p], 



(9) 



with 



A ± [p] = A (pQ [p]±e ± d (p A (p [p]\ , p% 



£ T ±((p-(p 0 ) 



Making use of equation (7) applied to the binary probability 
distribution p^, we obtain 



£+£_ 
N 



(10) 



To calculate the above bound it suffices to find the 'distances' £+ of 
the channel from the boundary measured along the tangent line. 
For extremal channels, £+ = 0 and the bound is not useful. For non- 
extremal channels, the above construction will yield a finite F c \ 
provided that £+ > 0, that is, A^ can be decomposed into a mixture 
of channels lying on the tangent. Channels that have this additional 
property will be called 'cp-non- extremal'. They obey the standard 
precision scaling, and include all full- rank channels (that is, chan- 
nels lying in the interior of the SQC) 41 ' 42 . An explicit method for 
calculating £+ and hence the bound for a general quantum channel 
is described in the Methods. 

Channel extension. Even though the CS method is very general, 
there are interesting examples of '(p- extremal' decoherence models 
for which CS does not apply. In this case, one can resort to the more 
powerful but less intuitive CE method. 
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The action of the quantum channel can be described via its 
Kraus representation 47 : 

A<p[p] = I^(p)p^/(<P)> (11) 



with Kraus operators satisfying HjKj((py Kj((p) = t. Although this 
representation is not unique, different sets of linearly independent 
Kraus operators are related by unitary transformations 

K i ((p) = ^u ij ((p)K j ((p) y (12 ) 

j 

where u((p) is a unitary matrix possibly depending on (p. 
An equivalent definition of the QFI has been proposed 40 ' 42 

f q[ a <p [p]] = min ^(l^)), (13) 



where the minimization is performed over all (p differentiate puri- 
fications of the output state 

A 9 [p] = Ti E {\V 9 )C¥ 9 \}. (14) 

For pure input state, different purifications correspond to different 
equivalent Kraus representations of the channel, as in equation (12). 
For many quantum metrological models 40 , one can make an 
educated guess of a purification and derive excellent input- 
independent analytical bounds providing the correct asymptotic 
scaling of precision. Nevertheless, the method may be cumber- 
some especially when the channel description involves many Kraus 
operators. 

A simpler bound can be derived by exploiting the intuitive 
observation that allowing the channel to act in a trivial way on 
an extended space, can only improve the precision of estimation, 
that is, 

max F Q [A4p]]< max F Q [a^ ®l[p ext ]]. (15) 

P Pext 

This leads to an upper bound on F Ni which goes around the input 
state optimization, yielding 42 : 

^AT<4mn{iV|a^| + iV(iV-l)|^| 2 J (16) 

where ||-|| denotes the operator norm, the minimization is 
performed over all equivalent Kraus representations of A^, and 

<*K = ^((P^iitp)' Pk = ^H<P)U<P)> (17) 

i i 

where Kj((p) = d^K^cp). For any given (p= (p 0 , equation (16) involves 
only K x and its first derivatives at <p 0 . Moreover, the bound is insen- 
sitive to changing the Kraus representations with a (p-independent 
u. Therefore, it is enough to parameterize equivalent Kraus rep- 
resentations in equation (12) with a hermitian matrix h, which is 
generator of u((p) = exp[ - ih((p- <p 0 )]. This reduces the optimization 
problem equation (16) to a minimization over h. 

4 



As the SS of precision holds when J 7 ^ scales linearly with N, the 
bound equation (16) implies that a sufficient condition for SS is to 
find h for which fi^ = 0, or equivalently 42 : 

Jh lj K}K r ijK\K q . (18) 

i,j q 

Here we go a step further and show that in this case we can obtain a 
quantitative SS bound 



F N <4ATmin \\a^ , (19) 
h 

where the minimization runs over h satisfying equation (18), and 
can be formulated as a semi-definite program, as described in the 
Methods section. 

Moreover, it turns out that the bound resulting from find- 
ing the global minimum in equation (19) is at least as tight as the 
one derived using the CS method based on equation (7) (see 
Supplementary Methods) proving superiority of the CE over the 
CS method. 

Examples. All examples of channels presented below are of the form 
A<p[p] = A[U(ppU^], that is, a concatenation of a unitary rotation 
encoding the estimated parameter <p and an (p-independent deco- 
herence process. Consequently, K i ((p)=K i U ( p, where i<Q((p), Kj are 
the Kraus operators of A^ and A, respectively. The most relevant 
models in quantum- enhanced metrology belong to this class, but 
the methods presented may be applied to more general models 
as well. 

In what follows, we adopt the standard notation where 11 is the 
2x2 identity matrix and are the Pauli operators. We focus 

on two-level probe systems (qubits) sensing a phase shift modelled 
using a unitary U((p) = exp[(ia^(p)/2] — rotation of the Bloch ball 
around the z axis. In the case of atomic clocks' frequency calibration 
(p=bco-t with 5ft) being the detuning between the frequency of the 
atomic transition and the frequency of driving field, while t is the 
time of evolution. Even though in practice the parameter to be esti- 
mated is 8ft), we will consider that to be (p, to have a unified notation 
for atomic and optical models. In the case of a two-mode optical 
interferometer, U((p) is the operator acting on a single photon state, 
accounting for the accumulated relative phase shift <p between the 
two arms of the interferometer. 

We apply the methods to four decoherence processes encoun- 
tered in quantum-enhanced metrology: two-level atom depolariza- 
tion, dephasing, spontaneous emission and the photon loss inside 
the interferometer. These examples will provide us with the full pic- 
ture of the applicability of the methods discussed in the paper, as 
their cover all distinct cases from the point of view of the geometry 
ofSQC. 

Depolarization. Two-level atom depolarization describes an iso- 
tropic loss of coherence and may be visualized by a uniform Bloch 
ball shrinking (Fig. 4a) where 0<T] < 1 is the final Bloch ball radius. 
Its description involves four Kraus operators 

which makes it an example of a channel lying in the interior of the 
SQC. Using CS method, we infer the SS of precisi on and calculat e the 
'distances' from the boundary of the SQC £+ = yj{l-r])(l + 3r])/{2r]) 
(see Methods) resulting in the bound presented in Table 1. Applying 
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Figure 4 | Graphical representation of decoherence models. Two-level atom decoherence processes are illustrated with a corresponding shrinking of 
the Bloch ball (blue) for the cases of: (a) depolarization, (b) dephasing and (c) spontaneous emission. The estimated parameter (p represents the angle 
of rotation about the z axis, whereas r\ specifies the strength of decoherence and effectively the size of shrinkage in the x-y plane. In the example (d) of 
the lossy interferometer, cp corresponds to the additional phase acquired by photons travelling in the upper arm and r/ stands for the power transmission 
coefficient in each of the arms. 



Table 1 1 Precision bounds of the most relevant models in 
quantum-enhanced metrology. 



Channel 
considered 



Classical simulation 



Channel extension 



Depolarisation 

Dephasing 

Spontaneous 

emission 

Lossy 

interferometer 



V(1-77)(1 + 377)/477 2 VO-77)(1 + 2t7)/277 2 
^-r\ 2 lr\ ^-r\ 2 lr\ 
NA (1/2)^1-77/7? 

^-V/ri 



NA 



NA, not available. 

The bounds are derived using the two methods discussed in the paper. All the bounds are of 
the form Aq> N > (const/ V/V), where constant factors are given in the table. Classical 
simulation method does not provide bounds for spontaneous emission and lossy interfero- 
meter, as these channels are ^-extremal. For the dephasing model, it surprisingly yields 
an equally tight bound as the more powerful channel extension method. 



the CE method, it is possible to further improve the bound and the 
analytical result of the semi-definite optimization is also shown in 
Table 1. The only no n-zero elemen ts of the corresponding optimal 
h are h 03 = h 30 = y/( 1 - 77) ( 1 + 3t])/c and h u = -h 2 \ = -i(l + rfjlc 
where c = 2(l-7])(l + 27]). To our best knowledge, the bound has 
not been derived before. 

Dephasing. Dephasing is a decoherence model of two-level atoms 
subject to fluctuating external magnetic/laser fields. In graphical 
representation, it corresponds to shrinking of the Bloch ball in x, 
y directions with z direction intact (Fig. 4b). The canonical Kraus 
operators read 47 



(21) 



where 0<T] < 1 is the dephasing parameter. As it involves only two 
Kraus operators, it is not a full-rank channel and lies on the bound- 
ary of the SQC. It is, however, a non-extremal and more importantly 
(p-non - extrem al channel allowing for the CS construction with 

£+ = yJl-T] 2 It] (see Methods), yielding the bound given in Table 1. 

Most importantly, the above bound, is exactly the same as the 
one derived by Escher et al. 40 or by the CE method, where the mini- 
mum in equation (19) corresponds to h = CJj/^^/l — r/ 2 j. This proves 
that despite its simplicity the CS method may sometimes lead to 



bounds that are equally tight as the ones derived with much more 
advanced methods even for channels lying on the boundary of 
the SQC. 

Spontaneous emission. The well-known two-level atom spontane- 
ous emission model is described by the Kraus operators 



K 0 = 



1 0 

0 Vn. 



n-n 

0 



(22) 



with 0<t] < 1. Interestingly, for all t] this channel is extremal 48 , which 
means that the CS is not applicable. Nevertheless, the CE method of 
equation (19) can be easily employed. 

Substituting Kj((p) into equation (18), we find that the genera- 
tor h is fixed to h = * A o z - 77I). Consequently, there is no need 

for minimization over h and from equation (19) we automatically 
obtain an SS bound listed in Table 1, which to our knowledge has 
not been reported in the literature before. 

Lossy interferometer. To model 'loss' in an optical interferometer, 
a third orthogonal state at the output, — vacuum — resulting from a 
loss of a photon, needs to be introduced. The decoherence channel 
on a single probe (single photon) is a map from a two- to a three- 
dimensional system: 



(23) 



where t] is the power transmission coefficient for the light travel- 
ling through each of the two arms. Although the corresponding 
channel A^ is non-extremal, it is unfortunately (p-extremal and CS 
cannot be used. Still, CE method can be easily applied. The opti- 
mal bound (see Table 1) corresponds to h with non-zero elements 





fo 0 N 




r 0 0^ 






°1 


K 0 = 


0 0 




0 0 




0 






v 0 








0 

V 


0 



h 00 =-h n 



2(1-77) 



. This bound is asymptotically equally tight to 



the best bounds known in the literature 38-40 , proving again that CE 
method despite its simplicity is able to provide powerful results in a 
straightforward manner. 

In optical interferometric applications, it is common to use states 
of light with an unbounded number of photons, such as coherent 
or squeezed states 49 " 53 . At a first glance, it is not obvious that the 
model considered in the paper covers these situations. Formally 
speaking, the bound we have derived applies to input states with a 
total number of photons fixed to N. Notice, however, that in every 



NATURE COMMUNICATIONS | 3:1063 | DOI: 10.1038/ncomms2067 | www.nature.com/naturecommunications 

© 2012 Macmillan Publishers Limited. All rights reserved. 



ARTICLE 



NATURE COMMUNICATIONS | DPI: 10.1038/ncomms2067 



optical experiment what is measured in the end are photon numbers. 
If all phase reference beams are taken into account, we can regard 
quantum states of light as incoherent mixtures of states occupying 
different total photon number sectors 53 ' 54 . From the point of view 
of metrology, the QFI is then bounded from above by the weighted 
sum of QFIs for each of these sectors 30 . Within each sector we can 
easily apply our bound and, as the bound is linear in N, the effective 
bound will simply correspond to replacing N with N— the mean 
number of photons in all relevant beams used in the experiment — 
and as such can be automatically applied to experiments involving 
coherent and squeezed light. 

To demonstrate the practical relevance of the bound, it is 
instructive to compare its predictions with the actual quantum 
enhancement observed recently in the GEO600 gravitational wave 
detector 17 . Although a detailed theoretical analysis of the set- 
up from the perspective of the derived bounds is still underway, 
by reducing the essential features of the set-up to a simple Mach- 
Zehnder interferometer, we can already give a preliminary estimate 
on how far is the actual experiment from the optimal perform- 
ance. For the reported overall optical transmission rj = 0.62, the 
theor eticall y predicted maximal quantum enhancement amounts 
to a y/l-T] =0.62 factor re duct ion in estimation uncertainty com- 
pared with the classical 1/yfrjN limit. The reported experimentally 
observed reduction was a factor of 0.67, which is an indication that 
the experiment operates close to the fundamental quantum limit 
and any significant improvement is possible only if the optical loss 
is further reduced. 

Discussion 

Assessing the impact of decoherence on the maximum possible 
quantum enhancement is a crucial element in developing quantum 
techniques for metrological applications. The tools developed here 
allow for a direct calculation of bounds on the precision enhance- 
ment for arbitrary parameter estimation model where the decoher- 
ence process acts independently on each of the probes and may 
be represented with a finite number of Kraus operators. The CE 
method is more powerful and for the most relevant metrological 
models yields in the asymptotic limit of large number of probes the 
tightest bounds known in the literature. Yet the CS method, which 
may fail to provide equally tight bounds, provides an intuitive 
geometric insight into the absence of asymptotic HS in the pre- 
sence of decoherence. From the derived bounds, it is clear that if 
the HS was to be preserved for large number of probes N the level 
of decoherence would have to decrease with increasing N roughly 
as (l-7])~l/iV. This gives an estimate on the regime in which the 
quantum enhancement is quadratic as compared with the regime of 
constant factor improvement. This is clearly seen in Fig. 2 where this 
transition appears around N~l/(1 - rf) = 20 for rf = 0.95. As in most 
metrology applications N is larger by several order of magnitude, we 
expect that the SS scaling provides a reliable bound for the optimal 
estimation precision. 

An important question that was not addressed here is the satura- 
bility of the bounds. The CS method does not provide tight bounds 
in general. As for the CE method, we are not able to prove that the 
bounds we have derived are saturable in the asymptotic limit. We 
should stress however, that these are lower bounds on estimation 
uncertainties and an upper bound can always be found by choos- 
ing a particular estimation method. Therefore, if a certain strategy 
performs close to the derived bound, we can certify that it is not far 
from reaching the fundamental quantum limit, as illustrated by the 
results of the GEO600 experiment. 

Even though the minimization over purification method 40 yields 
in principle a tight bound, in practice there is no effective algorithm 
to find a global minimum, and the only way to convince oneself 
that one has achieved the global minimum is again to show that the 



theoretical lower bound on estimation uncertainty coincides with 
the performance of a particular estimation strategy. 

The methods discussed were focused on deriving useful bounds 
in the asymptotic regime of a large number of probes. Still, the 
bounds are valid (though weaker) for any value of N. We leave it 
for a future work to improve the bounds for finite N, which seems 
to be possible by using the CE method and relaxing the fig = 0 
constraint. 

Methods 

Proof of the CS bound. To simplify the reasoning, let us focus on the CSs that are 
constructed using discrete sets of quantum channels, {A J. Then, the considered 
channel's action of equation (4) can be rewritten as a ^-independent map acting 
on a larger input space 41 



p [P] = Zifyi A i tPl = ° [P ® °<p 



(24) 



where cfy = Ejp^,- | | represents a diagonal state in some basis, in which O 
is denned via 0[p] = 2 2 (A 2 ®E i )[p\ with E i [a]=(e i \o\e i ). To prove equation (7), 
we write the QFI for the N parallel use of the channel and bound it from above, 
that is 



F Q [A^[p]] = F Q [o^[p®^]] 

< F Q [p®<r| N ] = F Q [o$ N ] = N F Q [<7„] = N F cl 



(25) 



exploiting the monotonicity of the QFI under any parameter-independent 
quantum map 55 , here 

Calculation of the CS bound. The geometry of the space of channels and more 
specifically the (p-extremality are best viewed by using the Choi-Jamiolkowski 
isomorphism 56 ' 57 . Given a quantum channel A : £(>4 n ) h-> £(?4 ut ) acting from 
the space of density matrices on H in to density matrices on H out , one defines 

P A = A ® I[|*F)(*F|], where |*F) = ^^ Hin | i)® | i) is a maximally entangled state in 
H[ n ®H in . A is a physical channel (that is, trace preserving, completely positive 
map) if P A is a positive semi-definite operator, satisfying Tr?^ {P A } =1 . If {£,•},• are 
the Kraus operators of the A channel, we can write explicitly P A = E z - | K 2 )(K 2 |, 
wherein) = K i <g>t\ x F). 

We can now say that the channel A^ is (p- non-extremal, if it is possible to find 
a non-zero £, for which P A ^ ± ed^Pj^ > 0. See Supplementary Methods for an 

alternative formulation of the (p-non-extremality condition and its relation to the 
well-known non-extremality condition due to Choi 57 . In practice, if we want to 
make most out of the bound in equation (10), we need to find the maximum values 
of £+, for which P A ^ ± £+3<pP A(?) are still positive semi-definite operators. This is a 

simple eigenvalue problem and therefore the bound can be obtained immediately. 
Taking the dephasing model as an example, the Choi-Jamiolkowski 



isomorphism of the corresponding A^ channel P A ^ 
has a simple form: 



--if^KiU^KiU^ 



1 0 0 TJe^ 

0 0 0 0 
0 0 0 0 



Tje 



0 0 



1 



(26) 



It is easy to check that P A ^ + sd^P^ > 0 provided | e \< <Jl - r\ /J], hence using 
equation (10) we arrive at the bound given in Table 1. 

CE method as a semi-definite program. Here we show that the minimization 
problem in equation (19) can be formulated as a simple semi-definite program. 
Let the channel A^ be a map from a d r to a d 2 - dimensional Hilbert spaces 
with Kraus representation involving k linear independent Kraus operators 
(^x^i matrices). Consider the following block matrix: 



A 



M d2 o 

0 y/tt d2 

0 0 



Ft 



(27) 
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where 1 is a dxd identity matrix. Positive semi-definiteness of the matrix A is 
equivalent to the condition 

a k ^ i k i <tl dl . (2g) 

i 

Minimizing the operator norm || || is thus equivalent to minimizing t subject to 
A>0. Taking into account equation (18), the problem takes the form: 

min t, subiectto : A > 0, Y/z^idX- = iYicjX. 
h } q ( 2 9) 

Since K" 2 - = K 2 - - the matrix A is linear in h and the problem is thus a semi- 

definite program with the resulting minimal t being the minimal operator norm 
|| (Xg ||. For the purpose of this paper, we have implemented the program using the 
CVX package for Matlab 58 . 
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